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Abstract 

The Dlx5 and D/x6 genes encode homeodomain transcription factors essential for the proper development of limbs in 
mammalian species. However, the role of their teleost counterparts in fin development has received little attention. Here, 
we show that d/x5o is an early marker of apical ectodermal cells of the pectoral fin buds and of the median fin fold, but also 
of cleithrum precursor cells during pectoral girdle development. We propose that early median fin fold establishment results 
from the medial convergence of d/x5a-expressing cells at the lateral edges of the neural keel. Expression analysis also shows 
involvement of dlx5a during appendage skeletogenesis. Using morpholino-mediated knock down, we demonstrate that 
disrupted dlx5a/6a function results in pectoral fin agenesis associated with misexpression of bmp4, fgfSa, and! and msx 
genes. In contrast, the median fin fold presents defects in mesenchymal cell migration and actinotrichia formation, whereas 
the initial specification seems to occur normally. Our results demonstrate that the dlx5a/6a genes are essential for the 
induction of pectoral fin outgrowth, but are not required during median fin fold specification. The dlx5a/6a knock down also 
causes a failure of cleithrum formation associated with a drastic loss of runx2b and colWal expression. The data indicate 
distinct requirements for dlx5a/6a during median and pectoral fin development suggesting that initiation of unpaired and 
paired fin formation are not directed through the same molecular mechanisms. Our results refocus arguments on the 
mechanistic basis of paired appendage genesis during vertebrate evolution. 
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Introduction 

In vertebrates, appendages (limbs, wings and fins) show major 
structural and functional difierences, but they present remarkable 
similarities in their developmental mechanisms [1]. Genes known 
to play a critical role in the initiation, growth, and patterning of 
tetrapod limbs (e.g. Tbx, Hox, Fgf, Bmp and Shh) are expressed in 
comparable spatiotemporal domain in fins [2-5] and share similar 
functions [5-1 1]. Particularly, the tetrapod Dlx5/Dlx6 and teleost 
dlx5a/dlx6a genes are expressed in the apical ectodermal ridge 
(AER) of the developing limbs in mice [12-15] and in the pectoral 
fin fold (PFF) and median fin fold (MFF) giving rise respectively to 
paired and unpaired fins in zebrafish [16-19]. At early stage of fin 
morphogenesis, teleosts present an AER structurally homologous 
to the AER of tetrapods [1,20]. Later, the AER transitions into an 
elongated pectoral fin fold [20] . AER and fin fold structures have 
been demonstrated to be essential signaling centers during 
appendage specification and outgrowth in vertebrates [7,17,20- 
22]. Despite the fact that tetrapod studies have demonstrated the 
central role of Dlx5/ 6 genes in limb formation, the implication of 
dlx5a/ 6a genes in teleost fin development has been little analyzed 
beyond examination of expression patterns. 



Dlx genes code for an evolutionary conserved group of 
homeodomain transcription factors, related in sequence to the 
Drosophila distalless gene {dltj essential for distal appendage 
patterning in insects [23]. The Dlx genes have arisen from the 
ancestral dlt gene as a result of gene duplication events [24] . In 
tetrapods, the Dlx family consist of six genes organized into the 
Dlxl/2, Dlx3/4 and D/xJ/ff bigene clusters. In zebrafish, eight dlx 
genes have been reported among which six {dlxla/2a, dlx3h/dlx4b, 
dlx5al 6d) are arranged on chromosomes similarly to their tetrapod 
counterparts ^'^11^. Expression and functional analyses oiDlxS/ 
6 have demonstrated their key roles in the development of the 
nervous system, of craniofacial structures, of endochondral bones 
and of appendages [13,14,28-38]. Simultaneous inactivation of 
Dlx5 and Dlx6 in the mouse results in a limb phenotype similar to 
that observed in patients affected with split-hand split-foot 
malformation type I (SHFM-I) [14,34,37]. Altered limb develop- 
ment in Dlx5/ 6 null mice is associated with loss of Bmp 4, FgfS and 
Msx2 expression in the medial part of the AER [14,15]. The data 
indicate that the Dlx5/ 6 genes have a central role in vertebrate 
appendage formation. It was therefore of interest to examine the 
function of dlx5a/6a during the development of paired and 
impaired fins. 
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Here, we show that dlx5a/6a genes are required for the 
initiation of pectoral fin outgrowth and for median fin fold 
morphogenesis. Our results suggest distinct requirements for 
dlx5al 6a genes in paired and unpaired fin development. More- 
over, the analyses demonstrate that dlx5a/6a are implicated in 
cleithrum formation and suggest that dlx5a is involved in fin 
skeletogenesis. 

Materials and Methods 

Ethical statement 

All experiments were performed according to the guidelines of 
the Canadian Council on Animal Care and were approved by the 
University of Ottawa animal care committee (institutional licence 
#BL 235). AH efforts were made to minimize suffering; 
manipulations on adult animals were performed with the 
anaesthetic drug tricaine mesylate (ethyl 3-aminobenzoate meth- 
anesulfonate; Sigma-Aldrich, Oakville, ON, Canada). Embryos 
were killed with an overdose of the latter drug. 

Animal maintenance 

Zebrafish and their embryos were maintained at 28.5°C 
according to methods described in [39]. Wild-type adult zebrafish 
were kept and bred in circulating fish water at 28.5°C with a 
controlled 14-h light cycle. Embryos were collected at the one-cell- 
stage. A Narishige IM300 microinjector was used for microinjec- 
tion. Wild-type, controls, and injected embryos were raised at 
similar densities in embryo medium in a 28.5°C incubator. 
Embryos were treated with 0.0015% 1 -phenyl 2-thiourea (PTU) to 
inhibit melanogenesis. 

Morphollno-medlated knock down and rescue 
experiment constructs 

For morpholino-mediated knock down, we injected or co- 
injected in 1 cell-stage embryos, 1 nl of dlxSa and/or dlx6a 
morpholinos at a concentration of either 0.4 mM or 0.8 mM. The 
choice of morpholino concentrations has been determined 
performing injection or co-injection of dlx5a and/or dlxf>a in a 
range from 0.2 to 1.6 mM. Injection of 0.2 mM MOs did not lead 
to any ()l)\ i()us phenotype whereas injection of 1.6 mM MOs was 
lethal after a few hours post-injection. The 5 '-untranslated region 
of each dbc gene was used to design translation-blocking antisense 
MOs against each dh transcript. The following translation- 
blocking MOs were obtained from Gene Tools (LLC, Philomath, 
OR, USA): dlx5a MO S'-TCCTTCTGTCG.^ATACTCCAGT- 
CAT-3'; dlx6a MO 5'-TGGTCATCATCAAATTTTCTGC 
TTT-3'. 

The following splice-blocking MOs that target exon 2 excision 
were also designed to confirm the phenotype obtained using the 
translation-blocking morpholinos: dlx5a e2i2 MO 5'-TATTC- 
CAGG.4AATTGTGCGAACCTG-3'; dlx6a e2i2 MO 5'-AAAT- 
GAGTTCACATCTCACCTGCGT-3' (from Gene Tools, LLC). 
Although the efficacy of the dlx5a e2i2 morpholino was deemed to 
be insufficient by Talbot et al. [38], in our hands, injection of 
0.4 mM dh5a e2i2 MO caused a 57% decrease in the levels of 
dlx5a transcripts. We observed comparable fin phenotypes 
injecting either dlx translation- or splice-blocking morphohnos 
(data not shown). 

As controls, vvc iiijcctc-d water or 1.6 mM of Standard Control 
MO (Gene Tools) that targets a human beta-globin intron 
mutation that causes beta-thalessemia. (Gene Tools 5'- 
CCTCTTACCTCAGTTACAATTTATA-3'). 

To ensure specificity of the morpholinos, rescue of the resulting 
morphant phenotypes was performed by co-injecting the corre- 



sponding dbc5a/6a mRNAs mutagenized on the MO target site 
{dlx5a MO binding site T(ATG)ACTGGAGTATTCGACA- 
GAAGGA, mxXdlxSa sequence C(ATG)ACGGGTGTTTTTGA- 
TAGGAGGA; dlx6a MO binding site AAAGCAGAAAATTT- 
G(ATG)ATGACCA, amidlxtia sequence ATTGCAAATAAT 
ATG(ATG)ATGACCA). Vectors were linearized with Notl and 
mRNAs were synthesized using the SP6 mMessage mMachine kit 
(Ambion). 

In situ hybridization 

In situ hybridization on whole-mount embryos and cryostat 
sections were performed as previously described [40,41]. 

For whole-mount in situ hybridization, embryos from 12 to 
96 hpf (n>100 for each experimental groups and markers) were 
fixed in 4% paraformaldehyde (PFA, MUlipore) in phosphate 
buffer saline 1 x (PBS, Amresco) overnight at 4°C, dehydrated in 
methanol, and stored in 100% methanol at — 20°C. 

For in situ hybridization on cryostat sections, larvae (n = 1 6) were 
fixed in 4% PFA in PBS overnight at 4°C, washed in PBS and 
equilibrated in sucrose in PBS overnight at 4°C. The samples 
were then embedded and frozen in O.C.T compound (Sakura 
Finetek) and sectioned at 10 |im. After in situ hybridization, 
sections were mounted with coverslips and Aqua Poly/Mount 
(Polysciences) before imaging. 

The antisense mRNA probes were labeled with digoxygenin- 11- 
UTP (Roche) and synthesized from cDNA clones: dlx5a [42], dlx6a 
[29], bmp4 [43],j^a [44], msxB [45], msxC [45], andl [46], runx2b 
[43], collOal [47]. NBT/BCIP (Roche) was used as alkaline 
phosphatase substrate. 

Whole-mount TUNEL assay 

The whole-mount TUNEL assay was performed on 24 hpf 
embryos (control embryos n = 62, dlx5al dhcGa morphants n = 98) 
with the Apoptag peroxidase in situ apoptosis detection kit 
(Millipore) following the modifications described in [48]. 

BrdU assay 

Dechorionated embryos at 23 hpf were placed in embryo 
medium containing 10 mM BrdU (5-bromo-2'-deoxy-uridine, 
Roche) with 1% DMSO for 1 h at 28.5°C (control embryos 
n= 13, dlx5a/dlx6a morphants n = 28). Embryos were killed with 
an overdose of tricaine mesylate and fixed in 4% PFA in PBS 
overnight at 4°C, washed in PBDT (PBS DMSO 1% Tween 
().\%), dehydrated in methanol and stored in methanol 100% at 
— 20°C. Then, the samples were rehydrated in a graded methanol- 
PBDT series, treated with proteinase K (10 |a.g/ml) for 20 minutes, 
post-frxed in PFA 4% for 20 minutes and incubated in 2N HCl for 
1 h at room temperature. 

Proliferating cells were immunodetected in wholemount 
embryos using a 1:100 dilution of the mouse anti-BrdU 
monoclonal antibody (Sigma), a 1:200 dilution of secondary 
anti-mouse HRP-conjugated antibody (Jackson Immuno), and 
revealed with DAB chromogenic substrate (Abeam). 

Histological picrosirius red staining 

48 hpf embryos (n= 119) were frxed overnight in 4% PFA at 
4°C and washed in PBST. The embryos were incubated for 1 h in 
a 0.2% solution of sirius red (direct red 80, Sigma-aldrich) 
dissolved in saturated picric acid (Sigma-aldrich). The staining was 
followed by washes in distilled water. Once the water ran clear and 
without any red color, embryos were sequentially dehydrated into 
glycerol/PBS solutions and stored in 100% glycerol. 
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Results 

dlx5a/dlx6a genes are early markers of apical ectodermal 
cells In developing paired and unpaired fins 

Expression of the zebrafish dlx genes has been previously 
examined [16,19,38,49] including in the developing appendages. 
Even before the median fm fold becomes distinguishable, dlx5a 
transcripts are expressed in ectodermal cells underlying the 
periderm at the lateral edges of the neural keel at 15.5 hpf 
(Fig. 1 A, A'). From 15.5 hpf to 16 hpf, (//x5a-expressing cells 
follow a dynamic convergent movement toward the dorsal midline 
to form the presumptive median fm fold (MFF) (Fig. 1 A-D), the 
anterior expression limit corresponding to the MFF domain 
around the 8'*" somite [17]. Then, dlx5a expression is limited to 
MFF ectodermal ceUs at 24 hpf and 48 hpf (Fig. SI), and 
gradually decreases until 72 hpf when transcripts are hardly 
detectable (data not shown). Expression of dlx6a mirrors that of 
dk5a except that transcripts seem to be present at lower levels (Fig. 
S2 B, D), a difference that was observed by us and others 
throughout the embryo using a variety of probes for this gene 
[18,19,29]. 

At the pectoral level, dlx5a transcripts are first detected at 24 hpf 
in apical ectodermal cells of the presumptive pectoral fm bud 
(Fig. 2 A). At 36 hpf, dlx5a is highly expressed in the AER of the 
developing pectoral fin buds (Fig. 2 B). From 36 hpf to 48 hpf, the 
AER develops into the pectoral fin fold (PFF) in zebrafish embryos 
[20], During PFF establishment, dlx5a expression is maintained in 
apical ectodermal cells until 48 hpf (Fig. 2 C). Transcript levels 
progressively decrease in the PFF from 48 hpf to 72 hpf (Fig. 2 D). 
Similar observations were made for dlx6a, although transcript 
levels appear to be weaker (Fig. S2 A, C). 

Transcripts of dlx5a are also expressed in the developing 
cleithrum at the base of the pectoral fm bud which extends the 
AER/PFF ^//xjfl-positive domain ventro-lateraUy from 36 to 
72 hpf (Fig. 2 B-D). The cleithrum is one of the major bones of 
the pectoral girdle which supports the pectoral fins in bony fish. 



A 


P AER 

^cl 

36 hpf 


PPF 

24 hpf 


\ 
CI ^'^d 

^48 hpf 


D 
1 

72 


hpf 



Figure 2. Expression of dIxSa during zebrafish pectoral 
development. Whole mount in situ hybridization for dlxSa in the 
pectoral region of zebrafish embryos from 24 hpf to 72 hpf. During 
pectoral fin formation, dIxSa is expressed in apical ectodermal cells of 
the presumptive pectoral fin bud (PPF) at 24 hpf (A) and in the apical 
ectodermal ridge (AER) of the early pectoral fin bud (PFB) at 36 hpf (B). 
At 48 hpf dtxSa expression is detected in the pectoral fin fold (PFF) and 
weak expression is observed in endochondral cells (end) of the PFB. 
iVloreover, the transcripts are detected in the developing cleithrum (cl) 
from 36 hpf to 72 hpf (B-D). The dashed lines in B-C indicate the limit 
between the AER/PFF structures and cleithrum precursor cells. Note the 
absence of dlxSa expression in the PFF at 72 hpf (D). Scale bars 50 |im. 
doi:10.1371/journal.pone.0098505.g002 

Weak dlx5a expression is also detected in the endochondral disc of 
the pectoral fm bud at 48 hpf and 54 hpf (Fig. 2 C and data not 
shown). These observations show that dk5a, and to a lesser extent 
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Figure 1. Expression of dlxSa during early specification of median fin fold ectodermal cells. (A-C) In situ hybridization for dlxSa in 
zebrafish embryos from 15.5 hpf to 16 hpf: dorsal view of the posterior axis (A, C) and coronal section (A') at the level indicated by the dashed line in 
(A). At 15.5 hpf, dIxSa is expressed in ectodermal cells at the lateral edges (le) of the neural keel (nk) underlying the periderm (pe) (A, A') (the dashed 
line in A' delineates the neural keel). From 15.5 hpf to 16 hpf, ectodermal cells expressing dlx5a follow a dynamic convergent movement to form the 
presumptive median fin fold (MFF) at the dorsal midline of the embryo (B-C). (D) Schematic representation of the zebrafish dorsal cellular movement 
implicating dIxSa based on A-C. The convergent movement produced by the establishment of the neural rod (nr) (16-19 hpf) leads to the fusion of 
the two lateral edges at the midline into the presumptive MFF expressing dlx5a. ep, epidermis; nc, notochord; psm, presomitic mesoderm; tb, tail 
bud. Scale bars shown in C for A-C and in A' 50 |im. 
doi:1 0.1 371/journal.pone.0098505.g001 
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dlx6a, are early markers of ectodermal cells giving rise to PFF and 
MFF structures and of precursor cells of the developing cleithrum 
and pectoral endochondral disc. 

Knock down of dixSa/da leads to severe appendage 
defects 

To analyze the implication of dlx5a/ 6a genes in appendage 
development, we performed dlx5a and dlx6a knock down in 
zebrafish embryos using morphoHnos (MO). We performed micro- 
injections at the 1 cell-stage with one morpholino or co-injection of 
two dbc morpholinos at two different concentrations (0.4 mM and 
0.8 mM). Embryos injected with translation-blocking MOs against 
dhSa and/or dlx6a exhibited characteristic and reproducible 
moderate to severe phenot^pes compared to control embryos 
(Fig. 3 A). When observed at 48 hpf, a moderate phenotype is 
defined as presence of a "curved tail" and hypoplastic pectoral fins 
whereas the severe phenotype corresponds to a "curly tail" 
associated with agenesis of pectoral fins. Moreover, embryos with 
a severe phenotype arc generally smaller in size when compared to 
controls and display craniofacial malformations (Fig. 3 C). 

In single dk5a morphants (0.8 mM), 118/236 (50%) of the 
embryos exhibit the normal phenotype whereas 111/236 (47%) 
show moderate and 7/236 (3%) display severe phenot^'pes (Fig. 3 
B). Knock down of dlxda (0.8 mM) leads to mostly normal 
phenotype (269/286) and we obtain a low proportion of moderate 
(13/286) and severe (4/286) phenotypes. The co-injection oi dlx5a 
and rf&ffMOs (0.4 mM each) increases the rate of moderate (232/ 
817; 28%) and severe (227/817; 28%) phenotypes. When we 
increase the concentration of the injected MOs (0.8 mM each), the 
double knock down results mainly in severe phenotype (704/950) 
and low rate of moderate (137/950; 14%) and normal (109/950; 
12%) phenotypes. Co-injection of splice-blocking dlx5a/dlx6a MOs 
was performed (n = 627) and confirmed the results obtained using 
translation-blocking MOs. The phenotypes of single dlx morphants 
compared to double dlx5a/6a morphant embryos underlie the 
potentially redundant function of the Dlx paralogs in vertebrates 
[13,14,50-52]. Moreover, the results show that the phenotypes 
observed in double morphants are dose-dependent. When we 
doubled the concentration of the MOs from 0.4 mM to 0.8 mM 
each, the rate of severe phenotype also increases from 27.8% to 
74.1%. Based on the above observations, we performed subse- 
quent experiments by injecting dlx5a + dlx6a MOs at 0.8 mM and 
we considered the severe phenotype embryos as specimens in 
which the dlx5a/dlx6a knock down was more efficient. Thus, the 
dlx5a/dlx6a morphants analyzed in the study are embryos 
presenting a severe phenot^qDe. 

We performed dlx5a/dlx6a mRNA rescue experiments to test 
the specificity of the phenotypes obtained with the dlx5a and dbc6a 
MOs (Fig. 3 C). We mutagenized 5'-UTR sequences in the dlx 
mRNAs that were co-injected with the dlx5al 6a MOs to prevent 
MO binding. The results show that co-injection oi dlx5al 6a MOs 
and dlx5al 6a mRNAs (70 ng/ |il each) increases three times the 
proportion of normal phenotype and decrease almost twice the 
number of severe phenotype embryos compared to the dlx5a/6a 
double morphants co-injected with GFP mRNA (200 ng/ |a.l) 
(Fig. 3 C). The experiments suggest that these phenotypes are 
specific to the dlx5a/6a knock down. The phenotypic aspects, 
including craniofacial malformations, pectoral fin and MFF 
defects, not fully rescued in the embryos with mild and severe 
phenotype, can be explained by the aberrant ubiquitous dlx5a/ 6a 
overexpression which leads to mild developmental defects when 
injected alone in the embryo (data not shown). 



The dlx5a/6a genes are required for the induction of 
pectoral fin outgrowth and cleithrum differentiation 

To better understand the effects of dlx5a/ 6a knock down on 
zebrafish fin development, we analyzed the expression of different 
markers known to be involved in appendage specification and 
morphogenesis, including some that have been shown to be 
affected in Z)^,?/^ ' mouse embryos [14,15]. 

First, we examined expression of hmp4, msxB and msxC genes 
that are expressed in the presumptive pectoral fin liud of 24 hpf 
CT embryos, at a position where pectoral fin buds will appear a 
few hours later (Fig. 4 A, E, G). In 24 h^idlxSa/ 6a morphants, we 
observed that the expression of the analyzed genes is lost in the 
presumptive pectoral fin bud during pectoral fin specification 
(Fig. 4 A', £', G'). 

We then analyzed the expression of markers that are involved 
during pectoral fin morphogenesis. At 48 hpf, hmp4 transcripts are 
found in the pectoral fin fold (PFF), the underlying mesenchyme 
and in the whole pectoral fin mesoderm (Fig. 4 B), whereas _^/5a 
expression is only detected in the PFF (Fig. 4 C). At the same stage, 
andl is expressed in epitheUal cells of the pectoral fm bud, and in 
distal mesenchymal cells invading the fold [46] (Fig. 4 D). 
Moreover, Expression of msxB is limited to PFF cells and to the 
adjacent mesenchymal cells whereas msxC is only detected in 
mesenchymal cells underlying the PFF [45] (Fig. 4 F, H). In 
contrast, 48 hpf dlx5al 6a morphants exhibit a severe decrease or 
complete loss of expression of the analyzed genes in the pectoral 
region at (Fig. 4 F', H'). 

These results show that the dlx5a/6a knock down leads to 
severely impaired or aholish("d expression of genes implicated in 
early vertebrate appendage development, a phenotype character- 
ized by agenesis of pectoral fin buds. 

As previously mentioned, dlx5a is highly expressed in the 
developing cleithrum of wild-type embryos from 36 hpf to 72 hpf 
(Fig. 2 B-D). The cleithrum is a dermal bone [1] located at the 
base of the pectoral fin buds and is the first bone that mineralizes 
at the axial k-\c\ during early zebrafish development [1,53]. In 
36 hpf control embryos, dlx5a expression in the cleithrum is 
associated with runx2b expression (Fig. 5 A), an early/intermediate 
stage marker of osteoblast cUfferentiation. At 48 hpf, runx2b 
expression is maintained in the cleithrum and is also detected in 
the opercular and ceratobranchial-5 bones at the craniofacial level 
(Fig. 5 C). At the same stage, cells of the developing cleithrum 
highly express the collOal gene (Fig. 5 E), an intermediate/late 
stage marker of osteoblast differentiation in zebrafish [47,53,54]. 
Interestingly, we show that dlx5al 6a knock down leads to a drastic 
loss otmnx2b and collOal expression in the cleithrum at 36 hpf and 
48 hpf (Fig. 5 B-F), even in embryos presenting a moderate 
phenotype with hypoplastic pectoral fin buds (data not shown). At 
48 hpf, loss of runx2b expression in dlx5a/ 6a morphants is also 
observed in the opercular and ceratobranchial-5 bone precursors 
(Fig. 5 D), structures which develop from a rfZjcJfl-positive domains 
(Fig. S3) [38,49]. 

The present data demonstrate that dlx5a/ 6a genes are required, 
directiy or indirectiy, for osteoblast differentiation of the cleithrum 
in zebrafish. 

Knock down of dlx5a/6a impacts on median fin fold 
morphogenesis 

To study the effects of the dlx5a/ 6a knock down in unpaired fin 
development, we analyzed the expression of bmp4,fgf8a, andl, msxB 
and msxC genes in the median fin fold (MFF) of dlx5a/ 6a 
morphants at 16 hpf, 24 hpf (during MFF specification) and 
48 hpf (during MFF morphogenesis). At 16 hpf and 24 hpf, no 
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A Lateral view B 




Figure 3. Phenotypes obtained with different dlx5a/dlx6a morpKioiinos and mRNA treatments. (A) Phenotypes observed in dlx5a/6a 
morphant embryos. Lateral view of control (CT) and dlx5a/6a morphant embryos at 48 hpf. The dlx5a/6a gene knock down results in moderate 
"curved tail" and severe " curly tail" phenotypes compared to controls. The moderate and severe phenotypes are associated with hypoplasia and 
agenesis of pectoral fin bud respectively as shown in the pectoral region magnifications. Scale bar for all panels 100 |im. (B, C) The graphics show the 
percentages of normal (blue bars), moderate (green bars) and severe (red bars) phenotypes obtained at 48 hpf following injection of different dIxSa/ 
dlx6a MOs and dlx5a/dlx6a mRNAs. For each treatment, the number (n) of specimens analyzed is indicated and each experiment was performed at 
least 3 times. The B graph shows the following treatments: control embryos injected with H2O; control embryos injected with a control MO (1.6 mM); 
single morphants injected with either dIxSa or dlx6a MOs (0.8 mM); double morphants co-injected with dIxSa and dlx6a MOs at two different 
concentrations (0.4 mM or 0.8 mM each). The C graph shows rescue experiments: control embryos injected with GFP mRNA (200 ng/ |al); control 
embryos co-injected with dlx5a/6a morpholinos (0.8 mM each) and GFP mRNA (200 ng/ and embryos co-injected with dlx5a/6a morpholinos 
(0.8 mM each) and dlx5a/dlx6a mRNAs (70 ng/ |il each). 
doi:1 0.1 371/Journal.pone.0098505.g003 



notable changes in bmp4, fgfSa and msx expression in the median 
fin fold are observed in the morphants compared to controls 
(Table 1, data not shown). At 48 hpf, msxB is expressed in distal 
mesenchymal cells of the MFF whereas msxC shows a larger 
anterior expression domain (Fig. 6 A-B, black arrowheads). Both 
transcripts are also detected in the spinal cord. In dlx5al 6a 
morphant embryos, msxB and msxC expression is still present in the 
spinal cord, but only a few msx-positive cells are detected in 
mesenchymal cells in the distal part of the MFF (Fig. 6 A'-B', 
black arrowheads), whereas aberrant msxB expression is detected 
in MFF apical ectodermal cells (Fig. 6 A', red arrowhead). In 
contrast, the fgfSa expression initially observed in MFF ectodermal 
cells of controls at 48 hpf is not affected in the morphants (Table 1, 



data not shown). The latter results suggested developmental 
defects of MFF mesenchymal cells. 

The MFF of dk morphants presents a granular aspect which 
seems to be associated with absence of actinotrichia. The 
actinotrichia are non-calcified fibrils which develop in the early 
zebrafish fins. They act as a scaffold for the migration of 
mesenchymal cells and, later, for fin ray formation during 
zebrafish fin skeletogenesis. The actinotrichia defects in our 
morphants were further investigated by examining andl expression 
during posterior axis development. The andl gene is a member of 
the actinodin family which encodes proteins that are essential 
structural components of zebrafish fin rays and fin folds and 
required for actinotrichia formation [46] . In control embryos, andl 
expression is observed in epithelial cells of the fin fold at 24 hpf 
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Figure 4. Impaired expression of bmp4, fgfSa, andl and m5xgenes in the pectoral fin region of dlx5a/6a morphants. Whole mount in situ 
hybridization for bmp4 (A, A'-B, B'), fgSa (C, C), andl (D, D'), msxB (E, E'-F, F') and msxC (G, G'-H-H') at 24 and 48 hpf in dorsal views of control (A- 
H) and dlx5a/6a morphant (A'-H') embryos. At 24 hpf, bmp4, msxB and msxC genes are expressed in apical ectodermal cells of the presumptive 
pectoral fin bud (PPF) in control embryos (A, E, G). In dlx5a/6a morphants, bmp4 expression is lost or altered in the presumptive pectoral fin bud (A') 
and the msxB and msxC transcripts are hardly detectable (E', G'). In 48 hpf control embryos, bmp4 is expressed in the pectoral fin fold (PFF), the 
underlying mesenchyme and in mesodermal cells of the pectoral fin bud (PFB) (B). At the equivalent stage, fgfSa is detected in PFF ectodermal cells 
(C), and andl expression is observed in the distal mesenchyme and in epithelial cells of the PFB but not in the PFF (D). The msxB gene is expressed in 
the PFF and the underlying mesenchyme (F), and msxC is detected in the mesenchymal cells but not in the PFF (H). In contrast to what is observed in 
controls at 48 hpf, dlx5a/6a morphants show a marked decrease or loss of expression of the PFB markers associated with pectoral fin agenesis (B'-D', 
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F', H'). (I) Schematic representation of the pectoral fin bud at 48 hpf summarizing the expression of dIxSa and the analyzed PFB markers in their 
corresponding cellular types. Scale bars 50 |j,m. 
doi:1 0.1 371 /journal.pone.0098505.g004 



and in all MFF cells at 48 hpf (Fig. 7 A-B). As previously reported 
[46] , andl transcripts are not co-expressed with dlxSa in the apical 
ectodermal cells at 24 hpf, but rather in epithelial cells directiy 
adjacent to the apical ectodermal tissue (Fig. S4 A). In 24 hpf 
dlx5a/6a morphants, andl is detected in the MFF but transcripts 
show a shorter antero-posterior expression domain (Fig. 7 A', Fig. 
S4 B). At 48 hpf, transcripts are also observed in a less extended 
MFF domain (Fig. 7 B'). 

We then performed picrosirius red staining to analyze the 
cellular organization of the median fm fold and to reveal the 
presence of actinotrichia. In CT embryos at 48 hpf, actinotrichia 
are clearly visible (Fig. 7 C, the white lines indicate their 
orientations and lengths). In dh5a/ 6a morphants, actinotrichia 
do not form and the fin fold is smaller and granular. 

We used the zebrafish ET-37 enhancer-trap line [55], 
expressing GFP in migrating mesenchymal cells of the fin fold, 
to examine MFF cell migration following dlx5al 6a knock down. In 
54 hpf control embryos, migrating mesenchymal cells in the MFF 
are aligned and show filopodia orientated in die direction of 



CT dlx5a/6a MO 




Figure 5. Knock down of dlx5a/6a leads to a defect of cleithrum 
differentiation. Dorsal views of whole mount in situ hybridization for 
runx2b and collOal on control (A, C, E) and dlx5a/6a morphant (B, D, F) 
embryos at 36 hpf (A-B) and 48 hpf (C-F). In controls at 36 hpf, runx2b 
is expressed in precursors cells of the cleithrum (cl) (A) whereas 
expression is absent in the pectoral region of dlx5a/6a morphants (B, 
black asterisks). At 48 hpf, expression of runx2b and collOal is detected 
in differentiating osteoblasts of the cleithrum which supports the 
pectoral fin bud (C, E). The runx2b transcripts are also observed at the 
craniofacial level in the opercular (op) and ceratobranchial 5 (cb5) bone 
precursors (C). In contrast, dix5a/6a morphants show a drastic loss of 
runx2b and collOal expression in the pectoral region (D, F, black 
asterisks) and of runx2b expression at the craniofacial level (D, blue 
asterisks). Scale bar shown in F for all panels 100 |im. 
doi:1 0.1 371/journal.pone.0098505.g005 
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migration (Fig. 7 D, the white line indicates the direction of 
migration). In dlx5a/ 6a morphants, mesenchymal cells are 
disorganized and do not show filopodia (Fig. 7 D-D'). Moreover, 
the median fm fold phenotype observed in the double morphants 
is associated with increased cell death and decreased proliferation 
at 24 hpf (Fig. S5). 

dIxSa is expressed during fin skeleton formation 

Several studies in tetrapod species have demonstrated the role of 
D/xJ/ff genes in endochondral bone development [14,56,57]. The 
expression oi dlx5a/ 6a has also been reported in some craniofacial 
bones that develop by either endochondral or intramembranous 
ossification [38,49]. To further examine the potential implication 
of dlx5a in appendage skeletogenesis, we completed dlx5a 
expression analysis at later stages during fin bone formation. We 
performed in situ hybridization on sections of wild-type larvae from 
around 6 days (4.2 mm) to one month post fertilization (8.7 mm). 

Expression of dlx5a was first detected in the perichondrium of 
the hypurals in 6.2 mm larvae (data not shown) whereas no 
expression was observed in the forming radials. Slightly later, in 
6.6 mm larvae, dk5a expression is maintained in the parahypural 
and hypural perichondrium (Fig. 8 B, black arrowheads) and 
transcripts are also observed in maturing chondrocytes (Fig. 8 B, 
blue arrowheads). In the same specimen, dlx5a is expressed in cells 
surrounding the distal radials during radial segmentation (Fig. 8 C, 
black arrowheads), including in the zone of segmentation (ZS) 
(Fig. 8 C, green arrowheads), and in developing lepidotrichia of 
the anal fm (Fig. 8 C, orange arrowheads), but not in the proximal 
radials (Fig. 8 C). Expression of dlx5a in the ZS is still detected in 
7.3 mm larvae while radial segmentation is almost completed (Fig. 
S6 B, blue arrowheads) but not in 8.7 mm larvae after radial 
segmentation (Fig. S6 C, black asterisks), whereas expression is 
maintained but decreased in cells surrounding the distal radials 
(Fig. S6 C). Moreover, dlx5a transcripts are observed in maturing 
chondrocytes (Fig. 8 D, blue arrowheads) and in the flanking 
perichondrium of anal proximal radials in 7.3 mm larvae (Fig. 8 
D, black arrowheads). Later, dlx5a expression is detected in well- 
developed lepidotrichia of the dorsal fin in 8.7 mm fish (Fig. 8 E, 
orange arrowheads). The results suggest that dlx5a is implicated in 
the formation of skeletal components of the developing zebrafish 
fins. 

Discussion 

dlx5a/6a genes are essential for the initiation of pectoral 
fin development 

During zebrafish development, the three dlx bigene clusters, 
dlxla/2a, dlx3b/4b and dlx5a/6a, are expressed in fin primordia 
[16]. However, no fin phenotypes have been observed in dlxla/2a 
and dlx3b/4b morphants (E. Heude unpubhshed observations, 
personal communication from A. Fritz, [30,58]). In mouse, 
combinatorial mutations of several Dlx genes have shown that 
only Dlx5/ 6 and Dlx2/ 5 mutants display distal limb malformations 
[14,36,59]. Thus, the Dlx5/Dlx6 genes seem to play a critical role 
in the formation of appendages in vertebrates. 

In zebrafish, the dlxSa/ 6a genes are early markers of apical 
ectodermal cells of the presumptive and developing pectoral fin 
buds and median fin fold (Figs. 1, 2, Figs. SI, S2) suggesting then- 
role in early specification and morphogenesis of these structures. 
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Figure 6. Impaired median fin fold expression of msxgenes in dlx5a/6a morphants. Whole mount in situ hybridization for msxB (A, A') and 
msxC (B, B') in lateral views of the posterior axis of control (A-B) and dlx5a/6a morphant (A'-B') embryos at 48 hpf. In controls, msxB and msxC genes 
are expressed in the spinal cord (sc) and in mesenchymal cells of the median fin fold (MFF) (black arrowheads A-B). Slight msxB expression is also 
observed in MFF apical cells (red arrowhead A). In morphants, msx expression is limited to a few distal MFF mesenchymal cells (black arrowheads A'- 
B') and aberrant msxB expression is detected in the MFF apical cells (red arrowhead A'). The dashed lines in (A') underlie the undulating and larger 
notochord (nc) displayed in the morphants. (C) Schematic representation of the dorsal median fin at 48 hpf summarizing the expression of dIxSa and 
msx genes in their corresponding cellular types. Scale bar shown in B' for all panels 50 ^im. 
doi:10.1371/journal.pone.0098505.g006 



At the pectoral level, dlx5a is expressed during pectoral fin bud 
initiation at 24 hpf and the dlxSa/ 6a knock down is associated with 
agenesis of pectoral fin buds (Figs. 3, 4). The data indicate that 
dlx5a/6a genes are required for pectoral fin bud induction and 
outgrowth in zebrafish. We show that dlx5a expression in the 
pectoral fin fold progressively decreases from 36 hpf to 48 hpf to 
become hardly detectable at 72 hpf (Fig. 2). Interestingly, this 
period corresponds to the AER/PFF transition followed by PFF 
elongation [20]. It has been shown that the unilateral surgical 
ablation of the PFF induces an activation of dlx5a expression in the 
reformed AER, whereas dlx5a is no longer detected in the control 
PFF in zebrafish [20]. The data suggest that dlx5a is essential for 
the first step of pectoral fin specification and morphogenesis until 
AER/PFF transition. In mouse, the Dlx5/ 6 genes are expressed in 
the AER of limb buds from E9.5 to El 2.5 [13-15]. However, the 
genes are not expressed in the limb field during limb bud induction 
at E8.5 (E. Heude unpublished observations, [60]). The data 
indicate that Dlx5/ 6 are not involved in the early specification of 
limb buds in mouse. Indeed, in contrast to what is observed in 
dlx5a/6a morphants, the Dlx5/6 nuU mutant mice do not show 
limb bud agenesis, the Dlx5/6 inactivation only alters limb 
morphogenesis. 



The dlx5a/6a knock down leads to severe decreases or loss of 
expression of bmp4, fgfSa, andl and msx genes. The Dlx5/6 null 
mutant mice show a loss of Bmp4, FgfS and Msx2 expression in the 
medial part of the limb buds [14,15]. A recent study demonstrated 
that the Dlx5 protein binds to conserved sequences in the 
proximity of the Bmp2 and Bmp4 loci in vitro, suggesting a direct 
regulation oi Bmp genes [15]. Moreover, BMPs appear to act as 
signaling relays between Dlx and Msx genes during ectoderm- 
mesoderm communication for proper limb development in the 
mouse [15]. Dlx5 seems to have a key role in the initiation of an 
ectoderm-mesoderm dialog during tetrapod limb morphogenesis. 
Our results suggest that dlx5a/6a expression in ectodermal cells 
acts upstream of bmp4,fgf8a and msx genes in the pectoral fin bud 
primordia to induce pectoral fin outgrowth in zebrafish. The 
dlx5a/6a genes seem to be crucial to induce the ectoderm- 
mesoderm communication at the basis of appendage outgrowth in 
zebrafish. However, it appears that the role oi Dlx5/ 6 in the early 
development of paired appendages diverged between teleosts and 
tetrapods during vertebrate evolution. 



Table 1. Markers affected by the dlx5a/6a knock down during the early development of pectoral and median fins. 







Pectoral fin 




Median fin 






24 hpf 


48 hpf 


16 hpf 


24 hpf 


48 hpf 


bmp4 


A 


A 


U 


U 


X 


fgfSa 


X 


A 


X 


u 


U 


msxB 


A 


A 


u 


u 


A 


msxC 


A 


A 


u 


u 


A 


andl 


X 


A 


X 


A 


A 



A, affected; U, unaffected; X, not expressed. 
doi:l 0.1 371/journal.pone.0098505.t001 
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Figure 7. Defects in actinotrichia formation in tKie median fin 
fold of dlxSa/Sa morpKiants. Lateral views of the developing median 
fin fold (MFF) of control (A-D) and dlx5a/6a morphant (A'-D') embryos. 
Whole mount in situ hybridization for andl at 24 (A, A') and 48 hpf (B, 
B'). In control embryos, andl is expressed in epithelial cells of the fin 
fold at 24 hpf and in all MFF cells at 48 hpf (A-B). In dlx5a/6a 
morphants, and! expression is detected in the MFF but transcripts show 
a less extended antero-posterior expression domain at 24 hpf (A'). At 
48 hpf, morphants display a thinner andl expression domain in the 
ventral and dorsal part of the MFF (B'). The picrosirius red histological 
staining in 48 hpf control embryos reveals the MFF cellular organisation 
with the presence of actinotrichia (C, the white lines indicate their 
orientation and length). In the morphants, the MFF is smaller and 
granular associated with absence of actinotrichia (C). The ET 37-GFP 
enhancer-trap line reveals the MFF mesenchymal cells displaying 
filopodia which migrate along the actinotrichia at 54 hpf (D, the white 
line indicates the direction of migration). The dlx5a/6a knock down in 
ET 37-GFP embryos leads to impaired MFF mesenchymal migration, the 
cells are disorganized and do not show filopodia. Scale bars in A-B 
50 nm, C 20 nm, D 10 ^m. 
doi:1 0.1 371 /journal.pone.0098505.g007 



Distinct requirement for dlx5a/6a gene function during 
paired and unpaired fin development 

It has been shown that the early dorsal MFF apical tissue 
consists of epidermal cells covered by periderm [17,61]. Our 
results show that, at the posterior level, dlx5a is first detected in 
ectodermal cells underlying the periderm at the lateral edges of the 
neural keel, with expressing cells later converging towards the 
midline to cover the posterior neural rod (Fig. 1) [62,63]. From 
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Figure 8. dIxSa expression during unpaired fin siceietogenesis. 

(A) Overview of the posterior axial skeleton of a one-month-old 
zebrafish. Endoskeletal fin supports are colored in blue and red squares 
indicate the structures analyzed in (B-E). (B-E) Whole mount in situ 
hybridization for dIxSa on 10 |im parasagittal frozen sections of late- 
stage zebrafish from 6.6 mm (16 days post-fertilization) to 8.7 mm (1 
month post-fertilization). In 6.6 mm larvae, dlxSa is expressed in the 
perichondrium of the hypurals 1 to 3 (hi -3) and of the parahypural (p) 
(B, black arrowheads) as well as in maturing chondrocytes of hi and h2 
(B, blue arrowheads). The transcripts are also detected in cells 
surrounding the distal radials (dr) during radial segmentation (C, black 
arrowheads), notably in the zone of segmentation (C, green arrow- 
heads), and in the developing lepidotrichia (Ip, C, orange arrowheads). 
Later, dIxSa expression is observed in the perichondrium (D, black 
arrowheads) and in maturing chondrocytes (blue arrowheads) of 
proximal radials (pr) of 7.3 mm larvae. In 8.7 mm larvae, dIxSa 
expression is maintained in the well-developed lepidotrichia (Ip) (E, 
orange arrowheads). Black-brown patches in B-C, E are melanophores. 
Scale bars B, E 20 nm, C-D 10 |im. 
doi:10.1371/journal.pone.0098505.g008 



15.5 hpf to 16 hpf, we show that ectodermal cells expressing dlx5a 
are laterally connected to the neural ectoderm. Thus, we propose 
that the presumptive dorsal MFF territory expressing dlx5a results 
from the medial fusion of the lateral edges of the neural keel 
during neural rod formation (Fig. 1 D). 

We show that the dlx5a/6a knock down leads to impaired 
median fin fold development associated with msx and andl 
misexpression in MFF mesenchymal cells, increased apoptosis 
and decreased proliferation in the MFF at 24 hpf (Figs. 6, 7 A, A'— 
B, B', Figs. S4, S5). Moreover, similar to what was observed in 
andl I and2 double morphants [46], the median fm fold anomahes 
of dlx5a/ 6a morphants are characterized by defective mesenchy- 
mal migration and absence of actinotrichia (Fig. 7 C, C'-D, D'). 
As already mentioned, andl is expressed in cells adjacent to the 
(//.»:,5a-expressing apical ectodermal tissue in the MFF at 24 hpf 
Our results suggest that the dlx5a/6a genes may be required to 
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maintain proper andl expression leading to actinotrichia formation 
during early MFF development. 

The pectoral fin bud and median fin fold share similar 
embryonic components including apical ectodermal tissue and 
mesenchymal cells migrating along the actinotrichia to invade the 
folds (Fig. 4 I, Fig 6. C). Both structures express a similar set of 
genes in the corresponding embryonic tissues. Indeed, early dlx5a 
expression is detected in apical ectodermal cells of the presumptive 
MFF at 15.5 hpf and in the PFF at 24 hpf suggesting a similar role 
for dlx5a in the specification of both structures. However, whereas 
the dixjri/ 6a knock down leads to the early loss of expression of 
markers associated with fin bud agenesis at the pectoral level, we 
did not observe obvious differences in the expression of bmp4,fgf8a 
and msx genes in the MFF of the double morphants compared to 
controls at 16 hpf and 24 hpf (Table 1). Moreover, the median fin 
fold begins to develop in dlx5a/6a morphants and later shows 
developmental defects at 48 hpf In contrast to what we described 
at the pectoral level, dlx5a/6a genes seem to be required for proper 
morphogenesis of the MFF but not for its specification. 

In pectoral and median fin buds, dlx5a expression is associated, 
directly or indirectly, with a('ti\'ation oi hmp4,JgfSa, msx and andl 
genes. Fgf8 and Dlx proteins are also detected in the MFF apical 
tissue of sharks suggesting that ancestral molecular mechanisms 
implicated in unpaired fin development might have been 
established during early vertebrate evolution [22]. It is now well 
accepted that paired appendages evolved after unpaired append- 
ages [22,64,65]. Despite their different embryonic origins [66,67], 
it has been suggested that paired and unpaired fms use a common 
suite of developmental mechanisms, a hypothesis mainly based on 
expression analyses [17,22,68,69]. The latter studies support that 
ancestral mechanisms of median fin development have been co- 
opted for the development of paired appendages. In contrast, our 
observations show that dlx5al 6a knock down has a different 
impact on zebrafish pectoral and median fin development 
(Table 1). The early communi[:ation between MFF apical 
ectodermal cells and underlying structures seems to take place in 
dlx5a/6a morphants. However, dlx5a/6a expression in apical 
ectodermal cells may be required for the specification of the 
pectoral fin bud. Therefore, although development of paired and 
unpaired fms may use similar molecular mechanisms, when it 
comes to those involving dlx5a/ 6a, differences in timing, expres- 
sion territory or usage of target genes may underlie the profound 
phenotypic difierences that we observed. Alternatively, paired and 
unpaired fins may use different rf/xJfl/fe-associated molecular 
mechanism for their development. The differential effects of 
zebrafish mutations affecting other genes on pectoral and median 
fin development also support difierences in mechanisms (for 
review, [70]). 

Cleithrum formation requires dlx5a/6a expression 

We show that dlx5a is highly expressed in the differentiating 
cleithrum from 36 hpf to 72 hpf in control embryos (Fig. 2 C— D). 
The cleithrum extends from the base of the pectoral fin and forms 
the posterior edge of the gill chamber. The cleithrum is an 
ancestral component of the pectoral girdle not homologous to any 
bones in mammals. It is present in all bony fish ancestors 
(Osteichthyes) and stem-group tetrapods but is absent in living 
tetrapods [71], except frogs [72]. Apparentiy, loss of the cleithrum 
during vertebrate evolution corresponds to the appearance of neck 
structures for head mobility. In zebrafish, the cleithrum is a dermal 
bone from mesodermal origin [1,73] and is the first bone to 
mineralize at the axial level during early zebrafish development. 
Expression of dlx5a in the cleithrum is associated with expression 
of runxSb at 36 hpf and runx2b and collOal at 48 hpf, early/ 



intermediate and intermediate/late markers of osteoblast differ- 
entiation respectively (Fig. 5 A, C, E) [47,53,54]. The expression 
analysis indicates that dlx5a is an early/intermediate marker of 
differentiating cleithrum osteoblasts. 

In dlx5a/ 6a morphants, expression of runx2h and collOal is lost 
suggesting impaired cleithrum development (Fig. 5 B, D, F). 
Absence of cleithrum is not an indirect effect of the loss of pectoral 
fins in dlx5a/6a morphants. Many studies reported experiments 
that led to the absence of pectoral fins while the cleithrum was still 
present, often without any signs of dysmorphologies [10,11,74— 
77]. Moreover, dlx5a/6a knock down also leads to a loss ot runx2b 
expression in the opercular and ceratobranchial bones (Fig. 5 C, 
D), skeletal structures which both develop in a dlx5a/6a-poskive 
context at the craniofacial level (Fig. S3) [38,49]. Our results show 
that dlx5a/dlx6a genes are required for cleithrum formation in 
zebrafish. 

dIxSa is involved in fin skeletogenesis 

To further study the potential implication of dlx5a during 
zebrafish fin skeletogenesis, we extended dlx5a expression analysis 
to later stages of unpaired fin bone formation. The fin skeleton is a 
mix of endochondral and dermal bones; the radial and hypural 
bones which support and articulate the fins rays originate from 
cartilaginous structures (Fig. 8 A, blue structures), whereas the 
lepidotrichia (fin rays) are dermal bones [65]. We show that dlx5a 
is expressed in the developing endochondral and dermal structures 
(Fig. 8, Fig. S6). During endochondral fin skeletogenesis, dlx5a is 
expressed in maturing chondrocytes and in the perichondrium of 
hypurals and proximal radials (Fig. 8 B, D). We never detected 
transcripts in proliferating or resting chondrocytes. The results are 
consistent ^\ith whut is obser\'ed during long bone formation in 
tetrapods and suggest a conserved Dlx5 implication in chondrocyte 
differentiation among vertebrate species. 

In teleosts, the radials arise from a common mesenchymal 
condensation which later segments into proximal and distal 
components. It has been suggested that fish radial segmentation 
corresponds to tetrapod joint interzone formation and that both 
processes share similar spatio-temporal expression of genes 
between zebrafish, chick and mouse [78]. Our results show that 
dlx5a expression is detected in cells surrounding the distal radials, 
notably in the zone of segmentation (ZS) during radial segmen- 
tation suggesting its role in the latter process (Fig. 8 C, Fig. S6 B). 
It has been shown that Dlx5/ 6 genes are early markers of the 
presumptive elbow joint in chick developing limbs [79]. The study 
reveals that co-expression of Dlx5 and Gdf5, a gene known to 
regulate joint formation, corresponds to the initiation of elbow 
joint formation. Interestingly, gdfS is also observed in the ZS in 
6.6 mm zebrafish larvae [80] at the same stage as dlx5a transcripts 
(Fig. 8 C). Our data further support the potential implication of 
Dlx5 in appendage segmentation/joint formation in vertebrates. In 
parallel, we observed that dlx5a is expressed in the developing 
lepidotrichia of median fins (Fig. 8 C, E). Altogether, the results 
associated with what is observed in the cleithrum (Fig. 5) indicate 
that dlx5a is involved in zebrafish osteoblast differentiation. 
However, the role of dlx5a in zebrafish skeletogenesis requires 
further investigations. 

Conclusion 

Our results demonstrate that dlx5a/6a genes are necessary for 
the specification and outgrowth of the zebrafish pectoral fin buds. 
However, the dlx5a/ 6a genes do not seem to carry out the same 
role during the development of pectoral and median fins 
suggesting differences in the molecular mechanisms controlling 
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the early development of paired and unpaired fins. The origin of 
paired fins during vertebrate evolution is still controversial 
[22,66,67,69]. Our results can refocus arguments and may open 
new evolutionary perspectives on the mechanistic basis of paired 
appendage genesis in vertebrate species. 

Supporting Information 

Figure SI Expression patterns otdlxSa during zebrafish 
median fin development. Whole mount in situ hybridization 
for dlxSa on lateral view of the posterior axis of 24 hpf (A) and 
48 hpf (B) zebrafish embryos and on 10 ^m parasagittal frozen 

sections (A', B') at the level indicated in (A, B). At 24 and 48 hpf, 
dlx5a expression is limited to apical ectodermal cells of the median 
fm fold (MFF). nc, notochord; nt, neural tube; sc, spinal cord. 
Scale bars 50 ^m. 
(PDF) 

Figure S2 Comparison of dlx5a and dlxSa expression in 
the developing zebrafish fins. Whole mount in situ hybrid- 
ization for dlx5a (A, B) and dlx6a (C, D) in the pectoral fin bud 
(PFB) (A, C) and in the median fm fold (MFF) (B, D) of 30 hpf 
embryos. Expression of dlxja is detected in apical ectodermal cells 
of both pectoral and median de\'eloping fins (A, B, black 
arrowheads). Expression of dlx6a mirrors dlx5a expression, 
however dlx6a transcripts seem to be present at lower level (C, 
D, black arrowheads). Scale bars 50 [Lm. 
(PDF) 

Figure S3 The opercle, the ceratobranchial-5 bone and 
the cleithrum all develop in a dlx5a-positive context. 

Dorsal and lateral views of whole mount in situ hybridization for 
runx2h (A, A') and dlx5a (B, B') in the anterior region of 48 hpf 
zebrafish embryos. Expression of mnx2b reveals the opercle (op), 
the ceratobranchial-5 bone (cb5) and the cleithrum (cl) (A, A'), 
structures which differentiate at early stage of zebrafish develop- 
ment. The dlx5a expression analysis at equivalent stage shows that 
the three bones develop in rf/xJa-positive domains, explaining the 
loss of mnx2b expression in dlxSal 6a morphants at the pectoral and 
craniofacial levels shown in Fig. 5 D (black and blue asterisks), ba, 
branchial arches. Scale bars 100 |im. 
(PDF) 

Figure S4 Expression of andl in the median fin fold of 
control and dlx5a/6a morphant embryos. Ventral view at 
the posterior level of whole mount in situ hybridization for andl in 
24 hpf controls (A) and dlx5a/ 6a morphants (B). The ventral view 
reveals that andl transcripts are not expressed in the apical 
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Figure S6 Expression of dlx5a in fin skeletal compo- 
nents of late-stage zebrafish larvae. (A) Overview of the 
posterior axial skeleton of a one-month-old zebrafish. Endoskeletal 
fin supports are colored blue and red square indicates the 

structures analyzed in (B-C). (B, C) Whole mount in situ 
hybridization for dlx5a on 10 nm parasagittal frozen sections of 
7.3 mm and 8.7 mm late-stage zebrafish. As seen in Fig. 8 C in 
6.6 mm larvae, dbc5a expression is still detected in cells 
surrounding the distal radials (dr) in 7.3 mm larvae, including in 
the zone of segmentation (ZS) (B, blue arrowheads) when the 
radial segmentation is almost completed. Later, after segmentation 
(8.7 mm), dlx5a expression is maintained but decreases in cells 
surrounding the distal radials and is no longer detected in the ZS 
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